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^ : ABSTRACT 



We study magnetohydrodynamic (MHD) standing shocks in inflowing plas- 



^-^ ' mas in a black hole magnetosphere. Fast and intermediate shock formation is 

explored in Schwarzschild and Kerr geometry to illustrate general relativistic ef- 
fects. We find that non-equatorial standing MHD shocks are physically possible, 
creating a very hot plasma region close to the event horizon. Shocked downstream 
plasmas can be heated or magnetized depending on the values of various magnetic 
field- aligned parameters. Then we may expect high-energy thermal/nonthermal 
emissions from the shocked region. We present the properties of non-equatorial 
MHD shocks and discuss the shocked plasma region in the black hole magne- 
tosphere. We also investigate the effects of the poloidal magnetic field and the 
black hole spin on the properties of shocks, and show that both effects can modify 
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the distribution of the shock front and shock strength. We find for strong MHD 
shock formation that fast rotating magnetic fields are necessary. The physics 
of non-equatorial MHD shocks in the black hole magnetosphere could be very 
important when we are to construct the central engine model of various astro- 
physical phenomena. 

Subject headings: accretion, accretion disks — black hole physics — MHD — 
relativity — shock waves — plasmas — 



1. Introduction 

We consider a theoretical implication of a possible link between the conditions of mag- 
netohydrodynamic (MHD) shocks and the resulting shocked hot and/or strongly magnetized 
plasma region very close to the black hole event horizon. In a series of our previous investi- 
gations in the context of general relativity (see, Takahashi et al. 2002, hereafter TRFT02; 
Rilett 2003, hereafter R03; Fukumura 2005; Takahashi et al. 2006, hereafter Paper I), it 
has been shown that shock formation in MHD plasmas inflowing onto a black hole can be a 
physically plausible mechanism for creating very hot (i.e., Tj > 10^^ K for ions and Tg > 10^ 
K for electrons) or strongly magnetized plasma regions, which possibly could be associated 
with subsequent thermal/nonthermal high energy emission. These previous studies are the 
extension of the earlier works on hydrodynamical (HD) shock formation in black hole accre- 
tion flows (e.g., Chakrabarti 1990; Lu et al. 1997; Lu & Yuan 1998; Fukumura & Tsuruta 
2004, hereafter FT04). These authors point out that black hole rotation is also important 
for determining the shock strength. The slow magnetosonic shock formation in relativistic 
ingoing flows near a black hole was first explored by TRFT0 2 . By extending the work on 



shock formation in relativistic winds by lAppl &: Camenzindl (Il988l ). these authors solved 



general relativistic jump conditions in accreting plasmas. It was suggested that relativistic 
slow/fast MHD shocks can significantly heat up the postshock plasma (generation of hot re- 
gions), which may produce a sufficient amount of high energy radiation. R03 systematically 
examined the types of the preshock plasma mainly for slow MHD shocks. Fast MHD shock 
formation in inflowing plasmas was first studied in Paper I where a non-rotating black hole 
was primarily considered. 

These studies were conducted mainly for the equatorial accretion fiows. Therefore, a 
natural next step is to extend these studies to two dimensional, non- equatorial flows. Also, 
non-equatorial shock-heated region would be attractive as a possible high energy radiation 
source above an accretion disk in the central engine of active galactic nuclei (AGNs). There- 
fore, we investigate in this paper the polar angle dependence of various shock properties 
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(e.g., shock radius, compression ratio, number density, magnetization, entropy generation 
and so on). 

The goal of our current study is to suggest where in the parameter space a MHD shock 
can possibly occur for specified black hole magnetosphere models. That is, for a specified 
plasma source, we hope that our model will be able to predict whether or not the MHD 
standing shock formation is possible. Also, if the shock does develop, then we hope that 
our model will be capable of suggesting the physical nature of the shocked plasma. Solid 
understanding of non-equatorial shock formation can be useful for comparing our theoretical 
implications with future observations. Furthermore, from observations of some Seyfert nuclei, 
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hole rotation as well, although in this paper we will confine our attention to the black hole 
rotation slower than that of the magnetic field lines of the magnetosphere. 

In a black hole magnetosphere, the magnetic field geometry should in principle be de- 
scribed by the trans-field force balance equation (the Grad-Shafranov equation), which de- 
scribes the mutual interaction between the frozen-in plasma and the surrounding global 
magnetic fields. As the plasma generates the toroidal/poloidal current distribution, the 
poloidal/toroidal magnetic fields are produced and the magnetosphere is constructed. Al- 
though solving this trans-field balance equation is an extremely difficult problem, there have 
been some attempts in the past to study the steady-state magnetospheric configuration 
around black holes (see, for instance, Mobarry fc Loveless 1986 for Schwarzschild geometry 
and Camenzind 1987 for Minkowski geometry). iNitta. Takahashi. fc Tomimatsul (jl99ll ) ex- 
tended these studies to Kerr geometry in order to analyti cally study a rotating black hole 
magn e tosphere, and found dipole-like fi elds (see Figure 3b in lNitta. Takahashi. fc Tomimatsu 
199ll ) . iTomimatsu fc Takahashi! (1200 ll ) solved a vacuum magnetosphere with a thin equato- 
rial disk and obtained a dipole-like geometry in the disk region and a uniform field ge ometr y 
along the rotational axis. A similar magnetospheric structure was discussed by ( 120021 ). 
However, for the accreting plasma very close to the horizon, it was found that the magnetic 
field lines are almost radial in the poloidal plane (see, e.g., Hirotani et al. 1992; Komissarov 
2005; Uzdensky 2005). Therefore, in our current paper we mostly adopt a conical geometry 
because we are interested in the accreting plasma very close to the event horizon. However, 
we will also explore the effect of different magnetic field geometries. Figure [1] schematically 
illustrates a magnetic field geometry in the poloidal plane adopted in this paper. We explore 
the possibility of MHD shock formation along the poloidal magnetic field lines, which may 
illuminate the underlying accretion disk, with possible application to the central engine of 
AGN in mind. Although the plasma fiows can originate from the accreting material in the 
equator, the exact plasma source can be diverse, and we do not specify them in this paper. 
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From the standpoint of recent MHD simulations, the past several years have seen the 
long-awaited advance in numerical general relativistic MHD simulations to study the geo- 
metrical nature of the black hole magnetosphere and the dynamical properties of the plasma 
emersed in the global magnetic fields. These dynamical simulations will help to offer the 
means to predict the bo undary cond i tions and parameters governing the black hole magne- 
tosphere. For example, iKoide et al.l (120001 ) explored the roles of the global magnetic fields 
coupled with the a ccreting plasma where the shock formation is found in the equator, while 
Hirose et al.l (120041 ) analyzed the structure of various magnetic fields around a rotating black 
hol e. The astrophysical processes that infl uence the spin evolution of black holes are studied 
bv iGammie. Shapiro. &: McKinnevI (120041 ). Through long time-evolved MHD simulations, 
McKinnev fc Gammiel J2004h found five main subr egions of the black hole magnetosphere in 
a quasi-steady state at the final simulation time. iDe Villiers et al.l (120051 ) recently showed 
that the unbound outflows can emerge self-consistently in the axial funnel region where 
the large-scale magnetic field spontaneously ar ises. The structure of large-scale Poynting- 
dominated jets is discussed by iMcKinneyl (120061 ). On the other hand, in order to gain further 
physical insights, here we emphasize that parallel more analytic, steady- state investigations 
should be equally important, and that is a major goal of our current studies. 

From X-ray observations of accretion-powered AGNs, particularly Seyfert galaxies, it 
has b ecome evident that accret i on disks play a cruc ial role in the central high energy activities 



[e.g. 



Nandra fc Pounds! Il994l : iPounds et al.l Il994] ) . Rapid instrumental progress of recent 
X-ray observatories (i.e., Chandra and XMM-Newton) allows us to further investigate the 
detailed dynamics of accreting plasmas in the immediate environment of supermassive black 
holes hosted in these AGNs. It has been widely accepted that the gravitational energy of the 
plasma is viscously dissipated into thermal energy in the course of accretion, producing the 
signature of thermal emission in the observed spectra in many Seyfert galaxies. Furthermore, 
the presence of the power-law X-ray continuur n component suggests the existence of magnetic 
fields in/around the accretion disks (see, e.g.. lHaardt fc Maraschilll99ll . for an AGN model). 
Therefore, in order to better understand the nature of the accreting magnetized plasma near 
a blac k hole, it is important to consider the role of both general relativity and magnetic fields. 
Meierl ( 120041 ) proposed that the inner part of the ingoing flows may enter a magnetically- 
dominated, magnetosphere- like phase in some Galactic black hole candidates (GBHCs) (e.g., 
GRS 1915+105) and perhaps some low-luminosity AGNs (e.g., NGC 6251) as well. His 
model is based on a strong, large-scale magnetic field structure around a black hole. These 
issues point to the importance of investigating the basic physics of a high energy source 
in the black hole magnetosphere. One motivation for our current paper is thus to provide 
a scenario, in terms of MHD shock formation in inflowing plasmas, that generates hot or 
magnetized plasma regions very close (within a few gravitational radii) to a black hole. 
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The structure of this paper is as follows. In §2 we briefly summarize our previous work 
on MHD shock formation in a rotating black hole magnetosphere and describe the trans- 
magnetosonic property of the plasma. We also show general relativistic adiabatic MHD shock 
conditions for the possible parameter space. The parameter dependence of the shock-related 
quantities is explored in §3, where we study the nature of our shock solutions and shock- 
included global solutions. We particularly focus on the non-equatorial MHD shock formation, 
by examining the polar angle dependence and place constraints on possible allowed shock 
regions in the parameter space. In §4 we discuss our results in terms of a black hole-plasma 
system coupled to a global magnetic field lines. Brief summary and concluding remarks are 
given in the last section, §5. 



2. Assumptions &; Basic Equations 

2.1. Black Hole Magnetosphere and MHD Accreting Flows 

Wc consider stationary and axisymmetric ideal MHD accretion flows in Kerr geometry. 
The background metric is written by the Boyer-Lindquist coordinates with the c — G — 1 
unit, 

ds^ = ( 1 — j dt^ + dtdcp 

-^^d<j>' - ^dr' - me^ , (1) 

where A = — 2mr + a^, E = cos^ 9, A = (r^ + a^Y — a^A sin^ 9, and m and a 

denote the mass and angular momentum per unit mass of the black hole, respectively. We 
take m = 1 throughout this paper. The black hole event horizon is th = ^ + Vl — o?. 

In the context of general relativistic ideal MHD, the basic equations governing plasmas 
consist of: (1) particle number conservation law: {nu").^a where n is the proper particle 
number density and is the plasma four-velocity, (2) equation of motion: T"^^ = where 
T"^ is the energy-momentum tensor for plasmas and (3) ideal MHD condition: u^Fap = 
where Fa/s is the electromagnetic field tensor. We also denote the poloidal plasma velocity 
as Up = —{u^'Ur + u^ug). The energy-momentum tensor T"'^ is given by 

T"^ = nfxu'^u'^ - Pg'^^ + ^ (^F^^F/^^ + ^g'^^F^^ , (2) 

where = F^j.F'^^ and ji — {p + P)/n is the relativistic enthalpy, P is the thermal gas 
pressure and p is the total energy density. 
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For a stationary and axially symmetric ideal MHD plasma, there exist five conserved 
quantities along a field line: total energy E and angular momentum L of the plasma, angular 
velocity of the rn agnetic field line Qp, particle flux per magnetic flux t] and entropy S (see 
Camenzindlll986l . R03). The total energy and angular momentum of the plasma are 



E 
L 



Brf, 



(3) 
(4) 



where flp = —Ftr/F^r = —Fte/F^e (see iBekenstein fc Oronlll978l ) and rj = nup/ Bp. 



In the Boyer-Lindquist coordinates, the toroidal component of the magnetic field seen 
by a distant observer is defined by B^p = (A/S)Fej., and the poloidal one is defined as 
= -[g'''^i*Frtf + g'^%*Fetf]Gt^ where = {l/2)y/^eaf3^5F^^ is the tensor dual to 
Fai3 and Gt = gtt + gt^^F- The dynamical timescale of accretion is assumed to be much 
shorter than diffusion timescale so that the adiabatic prescription is used for the infalling 
plasma as P = Kp^, where F is the adiabatic index, K is related to entropy S. po = mpfi is 
the rest-mass density [i.e., p = po + P/(T — 1)], and rr ip is the particle rn ass. We introduce 



the so called entropy-related mass-accretion rate (e.g.. IChakrabartilll990l ) as 



Ai = mprjK 



N 



(5) 



the polytropic index is given by F = 1 + Note that M. = Q for the cold flow limit 

(because oi K = Q). 



The poloidal equation is given by (see, e.g., iTakahashi et al.lll990l ) 

p\l + ul)=E'[{a-2M')p-k\ , 



(6) 



with / = -{G^ + GtL)/[p^{M'^-a)], a = gu + 2gt^VLF + g^^Vt],, k = {g^^ + 2gt^L + guL^)/pl, 
L = L/E, pI = gf^ - gttg<j,^ and G^ = gt^ + g^^Vtp = g^^{^F - t^), where u = -gt^/g^^ is 
the angular velocity of a zero angular momentum observer (ZAMO). The relativistic Alfven 
Mach number M is defined as 



M2 



A-KpnUp A-nprjUp 



52 

p 



Br, 



(7) 



The toroidal magnetic field B^ can be expressed as i?^ = —AnrjEpwf- 



Any physical MHD accreting plasma ejected from the plasma source with a small 
poloidal velocity onto a black hole must become fast magnetosonic at the event horizon, 
going through three magnetosonic critical points: a slow magnetosonic point (S: r = rs 
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where Up = usw), the Alfven point (A: r = ta where Up = uaw) and a fast magnetosonic 
point (F: r = rp where Up = Upw)- At each magnetosonic point the plasma speed is equal to 
each magnetosonic wave speed. Here, usw? ^a w and upw are the slow , the Alfven and the 



fast magnetosonic wave speed, respectively (see iTakahashi et al.lll990l . for their definitions). 
When the preshock flow is cold (i.e., P = or /i = mp), the ejected plasma is super-slow 
magnetosonic; the slow magnetosonic point will be absent in the cold flow solution. 

The function / becomes zero when the numerator of the expression of /, + GfL, 
is zero. There are one or two radii saisfying this condition and we will denote this radius 
(r = ta) the "Alfven radius". When the plasma poloidal speed becomes equal to the 
Alfven speed {up = uaw or = a) at this radius, this point is called the "Alfven point". 
Otherwise, it is termed the "anchor point" where the toroidal component of the magnetic 
field is zero (M^ 7^ «)• In the next section, let us apply MHD shock conditions to these 
accreting plasma flows. 



2.2. General Relativistic MHD Shock Conditions 



In TRFT02 and Paper I, we discussed in detail general relativistic Rankine-Hugoniot 
shock conditions in Kerr geometry that apply to accreting plasmas: particle number con- 
servation [nM^Jsh = 0, energy and angular momentum conservation [T"^]sh = and 
magnetic flux conversation [*-F"^]sh^a = 0, where ia is the unit vector normal to the shock 
front. The square brackets denote the difference between the values of a quantity on the two 
sides of the shock front. 

The poloidal magnetic field in our assumption is explicitly given by Bp{r) = C / (V AS r^) 
where C is a constant. To investigate the effects of the non-radial poloidal magnetic field, 
the poloidal field is parameterized by 5: 5 = for a purely conical field and 5 > for an 
outward diverging field. The conical magnetic field with 5 = is adopted in most of our 
computations unless otherwise stated. In the following, for simplicity we set the shock front 
to be normal to the magnetic field lines, = (0,1,0,0). These conservation laws lead to 
the following jump conditions: 

niu\ = n2ul 5 (8) 



( 



VAS r 
1 



Ml 
M| 



1 + 

^(^2 - 1) 

1 + 



1 



(9) 



- 8 - 



and the relation: 



ififif = {fiEY[{a-2M')f-k]-(^-^^yM^. (10) 

where fi = fi/mp, E = E/rrip and 17 = Airrimp/C. The subscripts "1" and "2" respec- 
tively denote the preshock and postshock quantities. The field-aligned conserved quantities 
{E,L,Qf,v) ci-re continuous across a MHD shock. Then, our entropy-related accretion rate 
can be expressed as: 



rup (It) \ 1 + 

Note that M. is continuous in shock-free fiow regions, but it is discontinuous across a shock 
{M.I < M.2) because of the entropy generation, where the adiabatic shock is related to 
extremely inefficient cooling processes. That is, most of the heat generated at the shock 
front is carried away (advected) with the accreting plasma. 



It should be noted from the shock condition expressed by equations ([9]) and ffTOj) that 
the plasma energy E is coupled with the parameter 7) as fjE, while the angular momentum 
L is always coupled with E a.s L = L/E through functions / and k. Also, the enthalpy fi in 
equation frequently appears with 17 as fir), which can be expressed in terms of parameters 
fjE, L and Qp through equation (fTOj) . When the downstream plasma becomes extremely 
hot (i.e., large fi2 3> 1), the coupled parameter, r)E, acts as a controlling parameter for 
determining the shock solutions. The equation ([9]) also shows that the loss of kinetic energy 
leads to the gain of the (toroidal) magnetic and thermal energy across the fast MHD shocks. 
Thus, in our computations, we will use the coupled parameter f)E. 

We now introduce some shock-related quantities that are useful for examining the prop- 
erties of MHD shocks. First, let us define the magnetization parameter a, which denotes the 
ratio of the Poynting fiux to the net mass-energy fiux of the accreting plasma seen by the 
ZAMO (see TRFT02). From the ZAMO's standpoint, a can be written as 



^ ^ B^g^^iytp - uj) ^ e-ah 



where h = 5'**(1 — ujL) and e = 1 — LQp- 

The local strength of a shock is measured by the ratio of the postshock to the preshock 
number density, ^2/^1, as 

A^^. (13) 
ni 
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Since any physical shocks that satisfy the second law of thermodynamics must be compress- 
ible, we require A > 1 for physically relevant shock formation. Then, the postshock plasma 
is heated up. In order to quantitatively evaluate shock heating, we furthermore define a 
dimensionless postshock plasma temperature 6 as 

ksT^l-J^ (14) 
rup 1 + N ' ^ ' 

assuming the equation of state for an ideal gas (i.e., P = nksT where ks is the Boltzmann 
constant and T is the plasma temperature). Notice that © <S 1 for the cold flow (Pi = 
and //i = mp). Depending on the temperature ©, the postshock flow, heated by the shock, 
may be capable of emitting high energy thermal radiation, although in practice it is also 
necessary to consider non-thermal radiation due to the magnetic field. 

To include the radiative effects on MHD flows, which means the breakdown of ideal 
MHD approximation, we should formulate the general relativistic version of radiative (non- 
ideal) trans-magnetosonic flow solutions, where E,L,r],flF,-M. are not conserved along the 
magnetic field line. However, the task is very complicated, and that is beyond the scope of 
our present paper. Therefore, we will not include the radiative effects on the postshock hot 
flows and consider ideal MHD trans-magnetosonic flows. 



2.3. Shock-Included Trans-magnetosonic Accretion Flows 

In the following we consider the MHD shock formation around a Schwarzschild and 
slowly-rotating Kerr black hole (i.e.. < ujh < where ujh = uj{rji), sec also Paper I). 
Although the black hole accretes the surrounding gas inward by its strong gravity, there 
are several mechanisms in general for decelerating the accreting gas: the centrifugal force 
due to plasma's angular momentum, magnetic tension and magnetic pressure. For accretion 
onto black holes a shock can be formed through these obstacles. Detailed studies of this 
possibihty will be a focus of much of this section. We restrict ourselves to cold preshock 
plasmas which are initially injected from a plasma source (e.g., accretion disk/torus). We 
impose that the cold plasma starts falling onto the black hole with small poloidal velocity 
(i.e.. Up < Maw) from this injection point (I) which is located somewhere between the outer 
light surface (L) and the Alfven point (A). The exact location of the injection point will be 
important when considering the boundary conditions of the accreting plasma that should be 
specifled by some accretion disk model, but that is beyond the scope of our present paper. 

We should point out the rotational energy per total energy, LQp, is one of the useful 
parameters for predicting the degree of magnetization of the shocked plasma flow. In Paper 
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I we discussed that there is a finite range in the value for Lflp in order for the multiple 
Alfven points to exist for the Qp > ujh case; that is, (LVLp).^^^ < LQp < {Lflp)^^,^. As the 
value of LQp becomes larger, the plasma becomes more magnetized because in the energy 
equation the magnetic term (second term) becomes more effective as LQp increases. 

Here, we briefly describe a shock-included trans-magnetosonic accretion flow. The 
poloidal equation ([6]) specifies the velocity of a plasma fiow for a given set of initial pa- 
rameters. When a MHD shock develops, both preshock and postshock fiow solutions must 
pass through appropriate critical points: that is, the fiow must pass through the fast magne- 
tosonic points twice {outer and inner fast magnetosonic points). As discussed earlier, for a 
given set of conserved quantities, there exist two Alfven radii : outer and inner Alfven radii 
for a slowly-rotating black hole case. We consider that an upstream preshock solution 
goes through the outer Alfven radius while the corresponding downstream postshock so- 
lution passes through the inner Alfven radius. The outer Alfven radius is identified as the 
Alfven point (A) while the inner Alfven radius corresponds to the anchor point for the 
upstream solution. The anchor point is located somewhere between the outer fast magne- 
tosonic point Fout and the fast MHD shock location. When a shock forms inside the anchor 
point, the magnetic field line is refracted away from the shock normal. Such a shock is 
called a fast MHD shock. In fast MHD shocks, the preshock plasma must be super-fast 
magnetosonic while the postshock plasma must be sub-fast magnetosonic. This postshock 
plasma must pass through another fast magnetosonic point (F) again before reaching the 
horizon (H) (i.e., I ^ A ^ Fout Fast Shock —>■ Fjn —>■ H). When a shock forms outside 
the anchor point, the magnetic field line is fiipped over across the shock normal. Such a 



shock is called an intermediate MHD shock (see, e.g.. lHadalll994l : iDe Sterck fc Poedtsll2000 



2OOII . for its physical significance). In intermediate MHD shocks, the preshock plasma must 
be super- Alfvenic while the postshock plasma must be sub-Alfvenic. This postshock plasma 
must pass through another Alfven point (A) and a fast magnetosonic point (F) subsequently 
again before reaching the horizon (H) (i.e., either I — Aout Intermediate Shock 

Fin ^ H or I ^ Aout Fout Intermediate Shock Ajn Fin H). When a shock 
occurs right at the anchor point, it is called a "switch-on shock" for which the magnetic 
field line of the preshock fiow is radial. The shock condition 1^ determines whether a MHD 
shock can form for a given field-aligned parameter sets. In such a global accretion solution 
with a MHD shock, the preshock cold plasma [i.e., A1i(rsh,Mf) = 0] is connected to the 
subsequent postshock hot plasma [i.e., Al2('"sh, ^1) > 0] through the shock formation where 
the Alfven Mach number jumps from Mi to M2 (see Paper I). In the next section we system- 
atically explore the parameter dependence of the resulting MHD shocks by varying certain 
primary quantities. 
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3. Numerical Results 

In this section, the shock-included trans-magnetosonic accretion solutions are solved 
in wide ranges of field-aligned fiow parameters, where slowly-rotating black hole cases are 
considered; that is, the magnetic field lines in the black hole magnetosphere rotate faster 
than the black hole itself. In this case, there are two Alfven radii on the r-M^ plane. 
The locations are specified by the two field-aligned parameters, namely Qp and L. These 
parameters Qp and L have the minimum and maximum values for allowing the Alfven radii. 
These values depend on the geometry of magnetic field lines. In our assumed conical field 
configuration, these values depend on the polar angle 9. We must change the values of flp 
and L in the acceptable range in order to search for physical solutions. AC88 and R03 both 
have suggested that F = 4/3 is not a good approximation even for large radial velocity. In 
many fast shock cases, R03 has found that F ~ 3/2 (i.e., iV ~ 2). Following their claim, we 
choose F = 3/2 (A^ = 2) in this paper. We have checked that this choice is not qualitatively 
too critical for our end results. 

As to the choice of field-aligned parameter sets for computations, there are additional 
three restrictions for selecting our parameters: regularity conditions at the outer and inner 
fast magnetosonic points and the shock condition connecting preshock and postshock trans- 
magnetosonic fiows. Three out of five conserved parameters are related by these restrictions. 
Consequently, for a given black hole spin a and the magnetic field line with polar angle ^sh, 
the shock location Vsh is determined when the five parameters are selected in the acceptable 
parameter range, and the corresponding shock quantities are obtained. If we modify one 
of the parameters [i.e., the field-aligned parameters {E, L,r],QF) + geometrical parameters 
(a, Tsh, Osh, S)], to find another global shock solution, we must change at least one more field- 
aligned parameter. Therefore, to clarify the ^sh-dependence of MHD shock solutions, we 
change one parameter at a time in such a way that the degeneracy in shock solutions is 
removed. This allows us to systematically examine each parameter-dependence of the shock 
solutions. As explained in section §2.2, the shock condition in equation (Q is conveniently 
characterized by fjE and LQp, and therefore in the rest of the computations we will basically 
use fjE and Lflp as control parameters. 



3.1. Non-equatorial MHD Shock Solutions 

We first show in Figure [2] the possible shock region in the poloidal plane (rsh,6'sh) for 
different {E, LQp) sets, where we can uniquely obtain a point on the poloidal plane {rgh, Osh) 
for a set of {E^mp). (a) Vtp = 0.1, (b) 0.14, and (c) 0.2 for the a = case and (d) 0.21, 
(e) 0.22 and (f) 0.25 for the a = 0.5 case. Along each curve in these figures, either E or 
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LQp is fixed while the other parameter (either LQp or E) is varied, as indicated by the 
numbers. For a more detailed discussion of the global nature of the trans-magnetosonic 
accretion solutions, we select some representative shock solutions. The labels (l)-(3) denote 
the selected models: the model (1) near-equatorial shock, the model (2) mid-latitude shock 
with the same energy as (1), and the model (3) near-equatorial shock with larger energy. 
See Tabled] for the characteristic parameters for these models. As mentioned in §2 the shock 
location rgh is constrained by the shock condition in equation ([9]) and the trans-magnetosonic 
property. 

Here, we should mention that, in the black hole magnetosphere with a given black hole 
spin a and magnetic field distribution 6, we can obtain the cubic parameter space {E,LQf 
and flp), and the dependence of various shock solutions [e.g., X{Osh)] on 17 is not focused 
here since we find that it is very weak compared with the other parameters. This is because 
fjE ~ constant once L and flp are specified. That is, larger (smaller) fj and smaller (larger) 
E yields roughly the same shock solutions for a fixed LQp (see Paper I for details). We 
thus show the possible shock parameter space by E instead of "fjE and use a certain specified 
constant ?7-value to plot the diagram. Figure [2] is a map of this cubic parameter space 
onto the (rgh, Osh) plane, corresponding to a slice of such a cubic parameter space by the 
Qp =constant planes. Thus, under the considered magnetosphere, physically valid MHD 
shocks are restricted to this cubic space. 

As E increases, the shock location tends to shift further out for a given specific angu- 
lar momentum L. This is because the upstream plasma with very large E generally means 
large kinetic energy, too, making it easier for the fiuid speed to exceed the magnetosonic 
speed. Similar results are obtained in the HD shock formation (see, Lu & Yuan 1998 for the 
equatorial shocks and FT 2004 for the non-equatorial shocks). Then, the shock can form 
sooner (i.e., further out). In order for accretion to effectively operate near the mid to high 
latitudes, the specific angular momentum L needs to decrease to reduce the centrifugal force. 
The accretion cannot take place otherwise. We have confirmed in separate calculations that 
larger Qp can allow the MHD shocks to form closer to the rotation axis with Ogh as small as 
~ 5°. 

Next, let us examine the parameter space that allows the shock formation. In Figure E] 
we display (fast and intermediate) MHD shock solutions for various energy E, specific angular 
momentum, L and polar angle ^sh for fixed angular velocity of the field line Qp- It is first 
noted that the allowed shock region is continuous but limited in the parameter space. For a 
given set of {E,LQp) the shock strength A generally becomes strong er as increases (i.e., 
the closer to the equator, the stronger the shocks). This is also seen in the non-equatorial 
HD shocks (see FT04). The obtained qualitative pattern appears to be independent of the 
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choice of Qp and the black hole spin a (compare the left and right columns in Figure [3]). 
We find that there exists the allowed region for the shock formation in the parameter space 
spanned by {E, Lflp) for a fixed flp and a given magnetosphere model. Outside the allowed 
region (shaded region), physically valid shock solutions disappear for the two reasons: (i) 
the shock condition by equation ([9]) is no longer met as E decreases (see, e.g.. Paper I) and 
(ii) the inner fast magnetosonic point disappears as E increases. That is, as you slice the 
entire parameter space along a constant LQp-cmve as in Figure [3l there exists a definite 
border line between the shock allowed and forbidden regions determined by either condition 
(i) or (ii) depending on the ii^-value. Additionally, there is a constraint on the parameter 
space in terms of the existence of the Alfven point. Because the parameter space allowed 
by the requirement of the Alfven point is wider than those by the conditions (i) and (ii), 
the allowed shock region shown in Figure [3] satisfies all the conditions by (i), (ii) and the 
Alfven point. 

In general we have confirmed that the minimum value of ^sh is determined by the above 
condition (i) while the maximum value of ^sh by the condition (ii). For accretion to take 
place in higher latitude regions (i.e., smaller 6s^), the angular velocity of the field line Qp 
is forced to increase. Because the range of Lflp-vaAvLe is restricted by the existence of the 
Alfven points on the shock-included trans-magnetosonic accretion solution (see Paper I for 
details), flp must increase as L decreases in higher latitude regions for a shock to occur, 
where L (x sin^ 6'sh. Thus, the allowed region in Figure [3] topologically shifts a.s Qp changes. 
That is, the shock is allowed more in the mid-high latitudes as fij? increases. In other words, 
high latitude shock formation can only be allowed by large Qp. In this case, the shock tends 
to become stronger with increasing Qp. This feature, due to the rotation of the magnetic 
field line, is a new aspect inherent to the MHD shocks which would not exist in the HD shock 
properties. We have obtained diagrams similar to Figure [3] for different 17-values and find 
that the 17-dependence of the shock solutions is very weak. As to the dependence on energy 
E, it is apparent that the shock strength A, for instance, is very sensitive to the change in E 
[see, e.g., the models (1) and (3) along the solid curve]. On the contrary, changing E under 
the fixed LQp value has only a weak effect on shifting the polar angle ^sh- For a given E, 
the shock can occur over a wide range of angle 6^^ [see, e.g., the models (1) and (2) along 
the dashed curve]. 

The topological nature of the parameter space for a = 0.5 is essentially the same as for 
the a = case (compare the left and right columns in Figure [3]), in the sense that larger 
Qp allows stronger shocks near the higher latitude regions. However, for rotating black 
holes the shock region does not reach as high latitudes as for non-rotating ones. We have 
also found that the black hole spin a generally amplifies the shock strength as in the HD 
shocks. However, the coupling among various parameters in the presence of the magnetic field 
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complicates the situation. Around a rotating black hole, it is seen that the allowed region in 
the polar direction appears to be narrower compared with the a = case (compare the left 
and right columns in Figure [3]) • This is because there is a tighter constraint on the value for 
LVLp (not for r]E) that allows shock formation. This is explained by considering the presence 
of the Alfven point on the accretion solution. Because the location of the Alfven point is 
determined by LVLp-, as a changes, the allowed rang e of LVLp is varied (na rrowed down) for 



the Alfven point to exist on the accretion solution (ITakahashi et al.l 119901 ). This black hole 



spin effect for the Alfven point location is seen in Figure (31 for the allowed shock location. 

We present other various shock properties for the postshock flows in Figure IH (a) 
shocked plasma temperature 0, (b) postshock entropy-related accretion rate M.2 and (c) 
postshock magnetization a2- The shock solutions marked by (l)-(3) in Figure H] correspond 
to the models (l)-(3) in Figure Temperature strongly depends on E. With increase 
of E, the temperature can rise significantly [e.g., compare the models (1) and (3)]. The M.2 
vs. ^sh behavior is roughly similar to the vs. 6'sh behavior, except that there is a peak 
angle in the mid-latitude. On the other hand, the downstream magnetization (72 acts quite 
differently from the other two. With small increase of E, 02 decreases significantly [compare 
the models (1) and (3)]. That is, and 02 show a clear anti-correlation, as expected from 
the conserved energy equation ([3]). We also see a trend that the equatorial shocked plasmas 
are generally more magnetized than the non-equatorial ones. In separate calculations, we 
have confirmed that increasing will amplify the downstream magnetization (J2 while 
it will reduce temperature 0. Furthermore, in Figure for a fixed E the downstream 
magnetization 02 becomes larger as the I/fi^?- value becomes larger. Therefore, this is a good 
indicator for estimating the degree of magnetization of the shocked plasma. 

In §2.2 we assumed conical magnetic fields, which would be reasonable around a black 
hole. However, the magnetospheric configuration other than the 5 = case may be possible, 
where 8 parameterizes the cross-section of the magnetic flux. Therefore, we will here explore 
the effects caused by some other magnetic field configurations. Figure [5] illustrates the 
dependence of shock strength A on 5. The pure split-monopole field geometry on the poloidal 
plane is realized with 5 = 0, while the outward divergent poloidal field is parameterized 
by (5 > 0. We find that the MHD shock tends to become stronger as b becomes larger. 
This is because for 5 > the accreting plasma is more accumulated (or concentrated) 
toward a smaller cross-sectional area with decreasing r due to the tapered magnetic field. 
With increasing 5, the accreting plasma in such a field geometry would also become slower. 
Furthermore, the outward magnetic pressure would be more effective in that case. A larger 
deviation from a purely radial field configuration tends to not only allow the stronger fast 
shock A, but also enhance the shocked plasma temperature 0. Note that the increase in 
5 shifts the inner fast magnetosonic point inward, while the outer one remains almost the 
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same (jTakahashill2002l ) . The valid shock location must lie between these outer and inner fast 
magnetosonic points. Since the location of the fast magnetosonic points shift as 6 changes, 
the allowed shock location also changes. 



3.2. Properties of the Shock-Included Trans-Magnetosonic Solutions 

In this section we consider the global, shock-included trans-magnetosonic solutions for 
the selected cases [i.e., models (l)-(3)]. We present in Figure [6] several representative solu- 
tions, where the models (l)-(3) are selected from Figure [2b. That is, the model (1) near- 
equatorial shock, the model (2) mid- latitude shock with the same energy as the model (1) 
and the model (3) near-equatorial shock with greater energy. These inflows all start from 
several Schwarzschild radii. 

The model (1) shows a relatively small downstream kinetic motion [i.e., 1^2 ('")| and 
1*2 (r)] whose magnitude almost remains the same until reaching the event horizon, meaning 
that the external forces are balanced with the gravity (see Figure [6^ and b). For the total 
energy to be conserved, a large fraction of the total energy then needs to be redistributed 
to somewhere else. In the model (1), the major part of the energy is redistributed to the 
magnetic field across the shock, and the postshock magnetization cr2(r) is relatively large 
((72 ~ 1.1 in Figure [Hli), while the shocked plasma temperature is not high (0 ~ 0.5 in 
Figure H^). Such an anti-correlation is expected from the shock properties shown in Figure |H 
Note that |-B(/,i(rsh)| < \B^2{rsh)\, as expected, for the fast MHD shock. 

The model (2) has the same energy as the model (1) except for a smaller angular 
momentum L at = 40°. Since the toroidal magnetic field near the event horizon scales 
as Bff, ~ B^^H oc sin^ 9h where the subscript H denotes the event horizon, smaller |-B0,2| 
is required for the model (2) as seen in Figure In this case, the magnetic component 
of energy is accordingly smaller, and the postshock magnetization (J2{r) behaves as shown 
in Figure Oi. To compensate for a smaller magnetic energy component, postshock kinetic 
energy becomes relatively larger (see Figure Ek). As a result, heating is not significant 
(G ~ 0.3 in Figure H^). We also find that the shock location rsh shifts outward and in higher 
latitude regions as in the HD case (see FT04). For example, as ^^sh changes from 81° in the 
model (1) to 40° in the model (2), the shock location shifts outward from Vgh = 2.3 to 2.7 
(see Tabled]). This is because the centrifugal force barrier becomes stronger as 6'sh decreases. 

As we discussed in the previous sub-section, in Figure |6] the value of LQp is larger in 
the model (1) than that in the model (2) for a fixed E. The downstream magnetization 
parameter (T2{r) in Figure Etl is thus greater in the model (1), meaning that the model (1) 
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shows the more magnetized downstream plasma. 

For the model (3) we choose a larger energy, £' = 7.2 to see the energy dependence of 
the shock solution. For comparison, the values for the angular momentum and the polar 
angle are roughly kept the same as in the model (1). Despite the larger total energy, the 
toroidal magnetic field strength at the event horizon \B^^h\ has a similar value to that in the 
model (1) because of the boundary condition at the event horizon. The energy at the shock 
front in this case is then redistributed more to the thermal and kinetic energies. In fact, 
we see that the downstream radial acceleration \du2{r) / dr\ is larger, and the downstream 
toroidal velocity Mf(r) is also larger, compared with the models (1) and (2). As expected 
from Figure H^, highly heated plasma flow (G > 2) is realized. 

Next, we will examine the energy distribution of the plasma to further understand the 
nature of the shocked plasma flows. In Figure [7] we show the components of energy of the 
plasma, namely fluid -Efluid [first term of the right-hand-side in equation ([3])] and magnetic 
-^magnetic [secoud term of the right-hand-side in equation ([3])] components, for models (1) 
and (2) in Figure In the course of accreti on the energy transp ort between the fluid 



(E'fluid) and the magnetic field (-^magnetic) occurs. iHirotani et al.l (119921 ) discussed the energy 
transport (or energy redistribution of the plasma) between the two components in shock-free 
accreting MHD plasmas. In a slowly- rotating black hole case (0 < ujh < ^f), the dominant 
component in E is fluid energy (i.e., -Efluid > -^magnetic), but the magnetic component can 
become comparable to or exceed the fluid component near the event horizon for lager Lflp- 
value. [see the model (1)]. Note that after the shock each component almost remains 
unchanged. In other words, the shock acts in such a way that this energy transport from 
the fluid to the magnetic field is even enhanced at the shock front. 



4. Discussion 



In order to explore a general trend for MHD shock formation for the vast amount of 
parameter sets, we sliced the entire parameter space which consists of all of the parameters 
{E,L,QF,V:''^sh,dsh',CL,6), by finding the allowed cubic parameter space consisting of three 
parameters {E, L,Qp), for given black hole magnetosphere models with the other two pa- 
rameters a and 6 fixed. Since it is very hard to explore all the parameter space spanned by 
all of the parameters together, our current approach greatly simplifies and yet help us better 
understand a comprehensive picture of the resulting shock solutions. Through our search 
for possible non-equatorial MHD shocks we find that the allowed shock region is constrained 
by various physical factors - especially the regularity conditions at the magnetosonic points 
and the shock conditions. Because the mathematical expressions for the regularity condi- 
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and fast magnetosonic points are not at all simple 
Takahashil |2002| ) . the topological appearance of the obtained shock regions (see, e.g., 



tions for the existence of the Alfven 
(e.g. 



Figure Ej) is complicated, but nevertheless we find it very useful. 

We find in general that stronger MHD shocks can form when the plasma energy is 
larger. With magnetosphere rotation the shock becomes stronger. The shocked plasma 
temperature has a clear anti-correlation with the shocked plasma magnetization. The more 
strongly magnetized plasma is formed for larger LQp. We also find that the energy transport 
between the fluid and the magnetic field can operate even more effectively across the shock 
front. This transition is internally caused by the transport between these components. Near 
the equatorial regions the shock is generally stronger. It is thus possible for the astrophysical 
accreting plasma that very powerful MHD shock formation takes place near the the equatorial 
event horizon. In the present calculations non-equatorial MHD shocks are possible up to high 
latitudes with 6sh as small as ~ 25°. We found in separate calculations that ^sh can be as 
small as ~ 5° when fii? is sufficiently large. If the poloidal magnetic field or resulting 
shock is dynamically unstable and switched on and off as the rotating plasma continues 
to fall onto the event horizon, some sort of quasi-periodic phenomena, perhaps with the 
dynamical frequency of n(rsh) may be expected (where fl = u'^/u^ is the angular frequency 
of the shocked fluid). On the other hand, particles (primarily the rmal electrons) would be 
accelerated at the shoc k front through first-order Fermi mechanism (jFermilll949l : lBarindll997l : 
Gieseler fc Jonesll2000l ) in the presence of the randomly-distributed, turbulent magnetic field. 
The shock front could then be partly re sponsible for generating the base of the relativistic 



winds/jets (e.g.. iQuataert &: Loebll2005l . for the stellar wind case 



When considering ingoi ng flows, we adopt rnostly a split-monopole poloidal field geom- 
etry near the horizon (e.g., iBlandford fc Znajekl 119771 ). to mimic the radial magnetic field 
lines very close to the horizon (see Figure 1). Based on the ideal MHD conditions, the 
ingoing plasma is assumed to be frozen-in to the magnetic field lines. The ingoing plasma 
we considered in this work needs to be injected from some plasma sources. Although we do 
not consider the exact origin of such sources, there are some possible sites for generating the 
ingoing plasma; e.g., the surface of the accretion disk and/or its corona. In this case, the 
magnetic field lines connecting to the black hole may come from the outer plasma sources 
at large radii. Although we consider the conical field lines near the horizon, the poloidal 
field line can also be bent toward the equator and be connected to the equatorial accretion 
disk at regions further away from the black hole. We expect that such regions are at least 
outside the potential well of the particle's effective potential. Note that in the description of 
the steady-state black hole magnetosphere, the global poloidal field coupling the black hole 
event horizon to the surface of a geometrically-thin disk has already been considered (e.g., 
see Nitta, Takahashi, & Tomimatsu 1991; Tomimatsu & Takahashi 1991 for the analytic 
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solutions, and Komissarov 2005; Uzdensky 2005 for the numerical results). The poloidal 
field configuration adopted in our model is based on these studies. 

In this paper, a typical value for the angular velocity of the field line Qp which allows 
MHD shock formation is about twice that of the Keplerian disk flKepi'f'ms) (where r^g is the 
the innermost radius). That is, if the black hole-disk connecting field lines are considered 
under the conventional Keplerian motion, adiabatic standing MHD shock formation is less 
likely to occur. Our results suggest that the shock may develop if the rotation of the magnetic 
field is super-Keplerian [i.e., flp > ^Kep('"ms)] due to some processes. Although we may 
specul ate on what this additional p rocess may be, that is beyond the scope of our current 
paper. iMcKinney &: Gammid (120041 ) compared near-equatorial stationary MHD infiows with 
their time-dependent numerical results between the innermost stable circular orbit and the 
event horizon. In their infiow solution, there is no MHD shock when the time-averaged 
numerical value of Qp is about 10% of the fixepl'^ms)- Their result is consistent with our 
condition for the MHD shock formation; that is, no shock solution exists for Qp ^ ^Kep('"ms)- 

We now show that the values of various primary parameters we chose are not arbitrary, 
but they are chosen with possible applications to the centr al engines of A GNs in mind. 



For instance, for a typical magnetic field strength of AGNs, iKrolik I (119991 ) estimates the 



magnetic field strength as -B ~ 4 x 10^ x G under the equipartition assumption where T5 
is the effecti v e tern perature of the disk in units of 10^ K. For radio-loud narrow-line quasars, 
Wang et al.l (I2OOII ) suggests that the field strength of at least 5 ~ lO'' G is required for 



magnetic-heated corona to operate. Following their implications, we took a likely value of 
Bp ~ 10^ G for a typical accretion-powered Seyfert nuclei {Mbh ~ 10'' M0) in our estimate. 
Some of the primary field-aligned parameters are L and 17. The physical value for the specific 
angular momentum L must lie in a certain range in such a way that the Alfven points do 
exist (see Paper I for details). Since we consider the trans-magnetosonic fiows, our choice 
of L should be reasonable, t] tells us the particle fiux per magnetic fiux. Using a commonly 
used value for the mass-accretion rate for Seyfert nuclei (from observational estimates). 
Mace ~ lO~^M0/year ~ 10^"^ gram/sec, we can estimate 17- value as ~ Mrc^/C"^ (in physical 
units) at the horizon. Note that we only take 1% of 10^^ gram/sec for "our" non-equatorial 
MHD plasma fiows. That is, Mpiasma ~ O.OlMacc ~ 10^^ gram/sec. This is because the 
plasma density of the non-equatorial fiows should be considerably lower than that of the 
equatorial fiows. Also, we have Bp ~ C/r^, rj = nru^c/Bp, and Mpiasma ~ r'^u'^nnip. By 
ehminating u"^, nip, r] from these equations, we get f) ~ r MpiasmaC^ I ~ 0.01, which is in the 
range of values adopted in our calculations. Note also that the estimated value of depends 
on the values of the black hole mass and the magnetic field strength. Depending on these 
values, ?7-value can be smaller or larger than 0.01 by orders of magnitude. 
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Also, in an attractive model of accretion-powered AGNs, the critical outer region of the 
disk where the inflowing plasma originates is also where the outflowing plasma is injected. 
For the outgoing flows, the Lorentz factor 7 of the observed relativistic jets/winds is generally 
E = Eoo ~ 7 ~ 10 (for microquasa r black hole systems) to 7 > 100 (AGNs and Gamma Ray 
Bursts (GRBs)) (e.g., lMeierll2003l ). That implies the presence of highly energetic plasmas 
being accelerated within the black hole magnetosphere. Because they are injected to form 
the jets/winds with such high observed energy, E^o, the ingoing plasmas also should initially 
possess the same order of magnitude of energy E ~ ii^oo at the time of their launch from the 
foot points on the accretion disk. Therefore, the original energy of the MHD flow can also 
be £' ~ 10 — 1OC0. A major part of the energy can be magnetic in this case near the plasma 
source on the disk. Our selected value for E [E > \) in our calculations is then justified. 
For completeness, we would like to mention that no plasma energetically bound to the black 
hole (i.e., E <1) appears to get shocked in our case. 

We are aware of the degeneracy of our solutions, i.e., the coexistence of shock-free and 
shock-included solutions for the same set of conserved quantities. To further investigate 
which steady-state solution between the two is physically accessible in nature, we will need 
to perform a dynamical analysis by considering a time-dependent perturbation in the flow. 
Such an analysis is beyond the scope of this paper. However, we may no te that some authors 
(e.g.. iFerrari et al.lll984l : iTrussoni et al.lll988l : iNobuta &: Hanawalll994l ) have already carried 



out such an analysis (no magnetic field in a flat space) and concluded that some shock- 
included solutions are indeed preferred (physically accessible by nature) to the corresponding 
shock-free solution. That finding probably justifies the motivations of our current work. 

Our current investigations carried out through the analytical methods can help gain 
deeper insight when compared with the results obtained by long dynamical MHD simu- 
lations. For instance, when the numerical experiments for black hole magnetospheres are 
studied under various initial conditions, the MHD shock front may be generated in some 
computational domain after a long calculation time to a quasi-steady state. Then, we can 
compare the field-aligned fiow parameters for MHD shock formation in our analytic stud- 
ies with the final values of the corresponding physical quantities obtained by the numerical 
simulations. 

Recently, long-time dynamical evolution of magnetized plasma accretion onto a black 
hole has been exten sively studied, and the global poloidal magnetic fields are obtained by 



some authors (e.g., iHirose et al.l |2004J : iMcKinney fc Gammid l2004j : iDe Villiers et al.l 12005 



^For a given LQ,p-Mdl\xe, the TyiS- value is nearly constant (see Paper I). In the present work we do find 
MHD shocks for £^ - 10 with 77 - IQ-^ and E ~ 100 with f] ~ 10"^. 
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McKinneyll2006l ) . The magnetic field confined in the funnel region is evolved to quasi-steady 



states at later times, while the magnetic field in the corona and disk fluctuates in magnitude 
and its direction. Our results may apply to some regions during such a quasi-steady-state in 
late simulation times. However, the accreting flows (whether equatorial or non-equatorial) 
are expected to be turbulent, and hence our analysis here, based on axisymmetry and sta- 
tionary assumptions, is not exactly applicable for explaining short timescale local turbulence 
in these regions. With this in mind, it is still important to compare our results with the fu- 
ture large-scale MHD simulations in order to better understand the physics in these complex 
regions. 



5. Summary and Concluding Remarks 

In Paper I we showed that the MHD shock formation can occur in the equator. In the 
present work we extended that work to the non-equatorial two-dimensional geometry. We 
systematically explored, for the first time, general relativistic MHD shock formation in such 
a geometry by employing the conserved quantities, and found the allowed shock regions. 
In summary we find that non-equatorial MHD shocks can form and could be a plausible 
candidate for generating a hot or strongly magnetized region over various latitudes in the 
non-equatorial plane. 

In various astrophysical objects, high energy activities are often associated with the 
magnetic fields along which the accreting plasma fiows. Our investigation of MHD stand- 
ing shock formation will be useful for our better understanding of complicated interactions 
between the plasma and the magnetic field in these situations. In the context of a strong 
radiation source required for the accretion-powered central engines of objects such as AGNs 
and GBHCs, the presence of non-equatorial shocks therefore can be attractive as a candidate 
for such a source, and therefore our current work may turn out to be very interesting for 
future observations of these objects also. 

Before closing, it may be emphasized that our search for the parameter space which 
allows MHD shock formation is carried out in a systematic manner in such a way that our 
choice of the parameter sets is not arbitrary. It is consistent, for instance, with the numbers 
relevant in application to, e.g., the environment of black hole magnetospheres around super- 
massive black holes in AGNs. Within this context, our results are quite general. However, 
it may be noted that in reality the MHD shock location in the actual astrophysical situation 
would probably trace a certain trajectory, rather than occupying the whole allowed shock 
region in the constrained parameter space, because the parameters allowing the shock forma- 
tion may be unique from case to case. Our hope is that our current studies, through analytic 
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and steady-state analysis, are very useful, because they can be complementary to those with 
time-dependent numerical simulations, in the sense that the latter can take advantage of, 
for instance, the constrained parameter space found with our current studies. 
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/ Ingoing Flows 




Equatorial Plane 



Fig. 1. — A schematic picture of standing MHD shock formation in black hole accretion 
flows where the poloidal magnetic field lines near the event horizon can be approximated to 
be radial (not drawn to scale). Note that non-equatorial ingoing MHD plasmas, originating 
from some plasma sources, are made possible in the presence of the poloidal magnetic fields. 
The Alfven point, outer and inner fast magnetosonic points (dotted curves) are labelled by 
A, Fout and Fin, respectively. 



Table 1. Characteristics of models (1), (2) and (3) with different {E,L) in Figure [2l 
a = 0,6 = 0,Qf = 0.14 and r) = 0.006 for all cases. 
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Fig. 2. — Acceptable parameter space for possible MHD shock formation on the rsh-^sh plane 
spanned by parameters {E,LQf). The filled circles (•) represent the fast MHD shocks while 
the open circles (o) denote the intermediate MHD shocks. We choose (a) Qp = 0.1, (b) 0.14 
and (c) 0.2 in the left column for a = and fj = 0.006, while (d) Qp = 0.21, (e) 0.22 and 
(f) 0.25 in the right column for a = 0.5 and r) = 0.005. Solid curves labelled by the thin 
numbers denote constant LOj^-curves, while dashed curves labelled by bold numbers denote 
constant ^-curves. The labels (l)-(3) in (b) denote the selected models. See Table [T]for the 
characteristic parameters for these models. 
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Fig. 3. — Shock compression ratio A as a function of polar angle 6*8^1 for various E and 
LQp corresponding to the solutions in Figure El The shaded regions represent the regions 
where no shock-included trans-magnetosonic solutions are possible (forbidden regions) for 
any {E,LQp) combination for a fixed {Qf,v)- The other notations are the same as in 
Figure [2l 
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Fig. 4. — Other shock quantities as a function of ^^sh corresponding to the solutions in 
Figure [2b: (a) Shocked plasma temperature 9, (b) postshock entropy-related accretion rate 
A^2, and (c) postshock magnetization parameter a2- The notations are the same as in 
Figure [2l 




Fig. 5. — The dependence of the shock strength A(rsh) on S for E = 6.1, L = 4.1 and a = 0. 
We set fj = 0.006 for 5 = 0, r} = 0.003 for 5 = 1 and fj = 0.001 for 6 = 2. The crosses (x) 
denote the anchor points. The intermediate MHD shocks are obtained to the right of the 
anchor points. 
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Fig. 6. — Shock-included trans-magnetosonic accretion solutions as a function of radius r: 
(a) the radial four-velocity u"^, (b) the toroidal four-velocity u'^, (c) the toroidal magnetic 
field B^, and (d) the magnetization parameter a. The MHD shock in the accretion solution 
{vertical downward arrow) is formed after the preshock solution passes through the outer 
fast magnetosonic point, Fout- The postshock solution then passes through the inner fast 
magnetosonic point, Fin, and falls onto the event horizon, H. The solutions labelled (l)-(3) 
correspond to the models (l)-(3) in Figure [2b. See Table [H for the model parameters. 
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Fig. 7. — Radial profile of the energy distribution of E for models (1) and (2) in Figure [2]d. 
(a) LVtp = 0.55 for the model (1) and (b) 0.23 for the model (2). See Table [T]for the model 
parameters in detail. 



